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ABSTRACT 

A multiple search and interpolator metn&cJ b** been developed in o'der to use the ’Telt*ut** model to map 
• fc ** ro *‘ nertia. soil moisture and cumulative de !y evaporation by p.art in§ from therm*.' and vis -bit data accused 
hr a» borne scanner. Interpolations are made n th n *nd between look-up tables gene r a: td prevoutiy by tne model 
wAich uses a reverse ten pe'ature simulation p xfJi 

Aher the necessary geometric. rad»o met- c and atmosphere corrections have been applied, the seiner data 
ate processed together with ground and me:eo“0*oc cal data to produce mo store and erase *r ion maps. Compari- 
son n . th available measurements indicates a tendency to underestimate soii moisture content; a smilar trend $?*hh 
tofae a'so verified for cumulative daily evaporation. Mapping of locally non-homogeneous areas encounters soma diffi- 
culties which are discussed. 


1. INTRODUCTION 

The Joint Research Centre Ispra is carrying out. within the "Tellus" project, an investigation in collaboration with la- 
Wacoiies frorr. the whole European Community on the apolication of"HCMM satellite thermal data for soil moisture and 
fceac budget evaluation »n zones of agricultural and environmental interest (1 J. In this framework, the so-called "Tell-us" 
Made? has been devetooed for bare and sc3rsefy vegetated areas, which uses separately the day and night temperatures Td 
and Tnat the same geographical point, rather than their difference Td*Tn. to determine thermal inertia, soil moisture and 
onwul.vt da.? y evaporation. 

The * - Te!l us* w model will not be presented and discussed here as this has already been done ( 1 j. The Daoer will rather 
diuAf the data processing methods used to handle the model in conjurction with the ground, meteorological and scanner 

acquired during a joint flight performed on September 13, 1977, to the north west of London, U.K. 

Tfenedaia are su.t-ble fo^ a detailed application of the model as a preparation to the use of HCMM data on a wider 
acafte. The p^rrose «s to reach t- ■» a point by point mapping of soil moisture and evaporation by start-ng from pixel by pixel 
i-*v and w«ct * ‘-mpe^ lures and albedo, processed together with ground and meteorological data varying on a 
fcefcSIvy fir’d basis or content ail over the flight strip. 

The process take into account many correcting procedures: correction of scan anc € effect on the spectrometrlc 
cy^t editor* of atmosphere effect, geometrical corrections on the night thermal data in order to achieve supposition 

Ar d jy uvbk thermal data. 

At a vciai;. moisture and evaporation maps are constructed and it will be possible to evaluate their reliability, rather 
hr reb:n»f u^?^es actually fgenera! trends and behaviour intra- and interfield) than for absolute values, due to the !ack of 
ground measurements at this stage of the research. OLviouslv, the data processing sesame is independent of the 
<**•»« «it-»formaiH** of the model and each further improvement of the latter should bring some corresponding improve- 
ment n the thematic maps produced. 


2. THE DATA AVAILABLE • BRIEF OUTLOOK OF THE JOINT FLIGHT EXPERIMENT 

Tfc*s experiment was carried out on September 13, 1977, on the Grendon Underwood Experimental Catchment. 
■•lA «.-%* of London. A data strip 8 km long and 1 .6 km wide was acquired from 1000 m altitude with a Daedalus OS 

scanner at night time minimum *5.43 am) and day time maximum (2.00 pm) sur'ace temperature. The ther- 
AilbioieltlvKpr^ is between 8 ana 14 the vis bte and near IR ba^ds are given in Tab‘e 1 . The ground resolution 
■s?Sn *f 9*m hr. Night and c'ay were clear and practically without clouds, with wind 0.2 m s at night passage and 1 .7 m s 
At day passage Irreasured a* 20 cm above ground level). 
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Simultaneously with the overflights, ground surface temperatures w e re measured with a PftT-5 radiometer with hand 
$*5*11 *5 #im and 8 core samples were taken in a ploughed field (bare toil} of Denchwoah clay at 1 5 cm depth* for mois* 
tore determination. 

Air temperature, humidity and pressure were measured at various altitudes with a tethered balloon, 40 km from the 
site, to be used for atmospheric effect determination. 


3. HANDLING OF THE "TELL-US" MODEL TO PRODUCE POINT BY POINT MAPS 

3.1 Search and interpolation procedures within the Te!!*us look-up table 

To apply the Tellus model, an algorithm was developed ( 1 ] which calculates, for a given set of input parameters, a 
graph in the night temperature and day temperature coordinates (Fig. la) made up of equi -thermal inertia and tcui-cela* 
five surface humidity curves crossing each other. Each couple of measured day and night temperatures results in a point 
contained within one of the irregular quadrilaterals of this lattice; the corresponding thermal inertia (THI) and relative 
surface humidity (RSH) values can then be determined by constructing the interpolated equi-THI and equi*R$H passing 
th ough the point. A second graph of equi-THI curves, in RSH and cumulative daily evaporation (ODE) coofd^ates 
I Fig. 1b), allows the calculation of CDE by starting from the already known RSH and the proper interpolated equi-THI 
corresponding to the THI value ca leu la ted in the first graoh. 

A set of routines has then been written, which perform the following sequence of actions for each (Td, Tn) couple of 
the thermal scanner data: 

localize the quadrilateral containing the point in Td. Tn coordinates (Fig. la); 

- construct the interpolated equi-THI and equi-RSH passing through the point; 

- determine THI and RSH values by linear interpolation between neighboring values; 

construct the useful segment of the interpolated curve in the second graph (Fig. 1b) by linear interpolation between 
the upper and lower bounding equi-THI; 

determine the CDE value by linear interpolation within this segment. 

The foresaid procedure leads in fact to search and interpolating routines within a complex table with two entries (RSH 
and THU and three values (Td, Tn, CDE) at each entry crossing. 

3.2 Various classes of input data - Interpolation between look-up tables 

The input data required by Tell us can be listed in the following way according to their various acquisition rates: 

Data valid for the who'e flight: weather data measured every hour: solar irradiance, windspeed, air temperature, air spe- 
cific humidity; initial surface temperature at midnight and sub-surface temperature at lower boundary depth (0.50 m) 

(soil) and upper boundary heinht (air). 

Data valid on a field by field basis: slope direction and slope dip; surface aerodynamic roughness; surface em'ss»vity . 

Data acquired pixel h\ pixel by the scanner: day and night surface temperatures; visible and near infra-red specuometric 
measurements, combined to produce albedo values. 

About the acquisition rate of the albedo values, it was verified in a preliminary study that, if the surface temperatures 
are to be used pixel by pixel, the same must be done for albedo which hos a great importance in the surface heat balance 
entering in the Tell us algorithm and which also vanes pixel by pixel within a single field. 

In the joint flight experiment the fields considered have very low or negligible slope dips so that this parameter, to- 
gether with slope direction was not considered. 

Surface aerodynamic roughness is another important parameter and the fields considered, where the Tell*us model is 
applicable or can at least be tested, have been grouped into three broad classes: ploughed bare soil with roughness esti- 
mated to 0.015 m; stubbles with roughness estimated to 0.020 m; burnt stubbles with roughness estimated to 0 010 m. 

Emissivity should be, at best, measured for each field surface type. Such measurements, however, are rot available for 
the flight over Grendon; it was estimated that the range of variation should be from 0 ^6 to 0 97 for the fields considered. 
As the model is hardly sensitive to this emissivity variation as separate input, a mean vaiue of 0.965 was assumed for the 
whole area. As a result of this analysis, the procedure described in section 3.1 had to be enlarged in the following way: 

• calculate with Tell us a senes of look-up tables by varying step by step albedo a$ in.»ut within suitable extrema; 

calculate this series of tables for as many input values of roughness as neces^ / .vit n the present case). 

The following sequence of actions is then performed for each (Td, Tn, A) pixel of mer data. A being albedo value: 

- retrieve the field in which the given pixel is localized and the roughness value at** ibuted to this field; 
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retrieve in the corresponding series the two tables with albedo input values on both sides of the A value; 

calculate THI and RSH from (Td, Tn> in each table as in section 3.1; 

determine the final T HI and ASH values by linear interpolation between the two values calculated in the two tables; 

* determine CD E as in section 3.1. 

The retrieval of the corresponding field identity and roughness value implies that a file has been constructed in which all 
the pixels in each of the fields considered land not in the other fields) are replaced by index values equal to the proper 
field number (field identity) and that a correspondence between field number and roughness value is available. 

The various search and interpolating processes mentioned before, will thus work, within and between tables generated 
by the Tell-us algorithm, on a pixel by pixel basis with four "channels": day temperature, night temperature, albeco and 
field index. 

Although the calculation of each of the tables by Te!t-u$ can be rather time consuming in prc portion to the accuracy 
requested (from a few seconds to a few minutes CPU time of IBM 370/165), the search and interpolating operations are 
very fast 119 min. for 312.000 pixels, i.e. 3.6 x 10‘ 3 sper pixel) and allow the use of Telhuswith pixel by pixel input. 


3.3 Mapping soil moisture and cumulative daily evaporation 

Cumulative daily evaporation is obtained directly by the process mentioned before and then maoped point by oo»nt. 
Soil moisture percentage in the 0 * 15 cm upper layer is obtained by using the linear correspondence betv«en thermal 
inertia and volumetric water content. Thermal inertia THI is indeed expressed as: 


THI*V-p e .h e 


where t)r heat capacity p. - (2 (1 — 6J + 4.2 0 ] 10* (J/m 5 /k) for not too high organic matter content [2 T . and :~e 
thermal conductivity h . - h Q + (h 3 & — h 0 ) ~ (W/m/k). 0 is the volumetric moisture content, s the pore vo:ume 
equal to 0.50 for the Oxford day which is almost present everywhere in the fields of interest. h 0 end are the appa- 
rent therms! conductivities for 0=0 and 0 = 0.5 respectively, with values 0.75 and 1.4 (W m V) |3j. It can be verified 
that in the practical range (0=0, 0 = 0.5] THI depends lineariy from 0 with fairly good approximation. 


4 . CORRECTION OF THE SCAN ANGLE EFFECT ON THE SPECTROMETRIC VALUES 


As the scan angle of the Daedalus instrument used for the flight is 77 c , the variation of the atmospheric thickness 
along the scan line causes a variation in the atmospheric effect from the centre of the line towards the edges, with a resjl* 
ting systematic effect in the dat? acquired, which must be corrected. Due to the difficulty to implement a cetermm sttc 
correction scheme which requires an accurate knowledge of the physical phenomena as a function of angle, the stat stical 
approach known as "along track averaging" was preferred and applied to thermal and visible data. 

The method is knov.n: all along the strip of acquired data, the mean values of the one pixel columns of data a'e cal* 
cufated from edge to edge of the strip, a polynomial fitting of the mean values along the scan line attempts then to des- 
cribe the general trend of the effect on the mean of the columns; if The effect is absent or negligee, the fitted c ur.e re- 
sults in a straight horizontal line. The difference between fitted curve and value at the central column element inac r poi^t) 
calculated at each column position, gives then the correction to be applied to the value of the element in th s pont on, in 
every scan line. 

In the present case, as a careful visualization of the data put in evidence a loss of contrast towards the edges, it was 
decided to correct also ths variance in a similar way. 

The correction curws exhibited a noticeable, although asymetrical, effect on mean and variance for day thermal and 
visible data and on the mean for night thermal data, while the effect on the variance is almost neg ig.ble in the latter case. 
See Figs. 2a to 2d and in particular the strong asymetry in Fig. 2c. 

The mean and variance corrections were applied at the same time for each channel with the following formula: 


/ - *o + b Ik. — x ) 


where x is the noncorrected value for element in scan line i and column i; x ‘ is :hc- corrected value for The same ele 
ment; x is the value of the fitted curve of the means at column j; x^ is the value of the mean at nadir (central column); 
b ( is the variance correction factor at column j calculated from the variance fitted curve. 
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*. CORRECTION OF THE ATMOSPHERIC EFFECT AT NADIR POINT 


The correction described above has the effect to give at each point die channel values it would have at nadir position 
horn die air plane, in the ce tre of the scan line. The data must next be corr ec te d from the effect due to the atmospheric 
layer at nadir. 

&1 Thermal data 


The intention was to correct them by using the air temperature humidity ami press ion data, acouired at various alti- 
tudes. in a model applying the radiative transfer equation to atmospheric layers* Unfortunately the data had to be collec- 
ted rather far from the site <40 km apart) and the calculated corrections resulted to be very questionable. 

Comparison with measurements made on the ground leads to the following situation: 

* night temperature average values are very close to the corresponding average scanner values acquired over the same 
zone: 4.0 * 0.6* C (PRT-5) instead of 4.1 ; C (Daedalus) for ploughed field and. similarly. 4.3 - 0.4 and 4.8 for a field 
with a mixture of burnt and unburnt stubbles, 4.9 = 0.5 and 5.65 for a tennis court. 5.2 = 03 and 5.8 for a road 
crossing; 

two series of ground measurements v/ere made for day temperatures with the same radiometer type (PRT-5) and at 
the same time approximately. The first series gave values systematically lower, from 2 to 5 and even 7 = , than the 
scanner values for seven investigated places; it is. however, conjecturable that this could be cue to a shortage in 
the power supplied by the portable batteries. The second series, where two places only were investigated, cave values 
in very good agreement with the scanner: 26.4 (PRT-5) instead of 26.8 (Daedalus) for ploughed field, 23.6 instead 
of 24.2 for grass of a sports field. Under such circumstances, where the atmospheric corrections for thermal data are 
difficult to be validly determined and likely to be rather small, it was decided not to take them into account, and to 
rely uniquely on the ground measurements, ft is therefore desirable that atmospheric sounding be performed in situ 
m future campaigns. 


S3 Visible and near infra red data 

As the sky irradiance was measured from the ground by an Exotech 100 radiometer in the region from 0.50 to 
1.1 pm. only the channels 4 to 10 of the Daedalus scanner were used with relative band widths given in Table 1 . It was 
assumed that the reflectance value for this wavelength range would be a reasonable approximation for the albedo value re- 
ferred to the whole reflective range. 

The atmospheric corrections were applied to each channel and the average reflectances of some fields were next cal- 
culated by using the sky irradiance value from the ground. The necessary calculations were obviously done to mtercalibrate 
the two radiometers and to take into account their band width differences. These reflectances were then comoarec to the 
reflectances of the same fields measured on the ground. Details on the atmospheric corrections and the comparison be- 
tween scanner and ground measured reflectances will be published later; a short outline will be sufficient for the present 
purpose. 

The transmittance r of each Daedalus channel which represents the atmospheric absorption factor, was estimated 
from the measurements made with the Exotech radiometer with a visibility of about 18 km and corrected for the scanner 
attitude by the factor r ;w /r x calculated with the Turner model (4|. 

The additive path radiance term L y for each channel was also obtained from the Turner model with 18 km visibility 
and corrected for the scanner altitude with data obtained from the LOWTRAN model |51- 

Both r and L factors are displayed in Table 1; as L. depends on the surface albedo, the values are given fo r the 
three soil surfaces investigated, they were calculated from ground measured reflectances. As it is seen, they are wry* close 
to each other, so that an average value over burnt and unburnt stubbles (more often found among the investigated fields), 
could be used everywhere to simplify the calculations. The reflectance is calculated by the following formula { 6 J: 


f =JL_ f ( 

9 irt 1 


ASl r > r, 


where r, and L P are respectively the transmittance and path radiance for channel i; ** • ,ver received by the ra- 

diometer in channel i. calculated through the calibration coefficients, from the data o> channel i; I is the sky irradiance 
for channel i measured from the ground; A is the area of the ground seen through the IFOV of the scanner; SI is the 
effective optical aperture of the radiometer. 
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As a global result, the comparison is the following, between ground measured ref tectances end scanner measured re* 
ftectances corrected from atmospheric effect: 

ploughed field - hare soil: 0.08S (ground) and 0.090 fscumer); 
burnt stubble field: 0.048 and 0.055; 

field e/dh mixture of burnt and non -burnt stubbles: 0068 and 0.082* 

They are in both cases average values over the corresponding field. The discrepancies range from 6 to 20 *- # which is st i 1 ! 
acceptable by considering that the Tell -us algorithm is hardly sensitive to such albedo variation. Variations above 40 or 
507, which occur from field to field and often within the same field, must however be taken into account. Decreasing :he 
albedo in input from 0.1 to 0.05 for instance has the effect to increase by about 10 : thermal inertia and soil moisture in 
output. 


& CORRECTION FOR THE SURFACE EMISSIVITY 

The surface emissivity e is used to pass from radiative temperature T.^ measured by the ground or the aircraft radio- 
meter, to the actual temperature T. with the help of the formula T = T,c~l derived from the Stefan-Bostzmann radia; on 
law. The e value varies with each field surface but as measurements are not available for the Grendon experiment e-*d esti- 
mates of € give a variation which is likely to be in the range of 0.96 to 0.97 for the fields of interest here, a constant ave- 
rage value of 0.965 was assumed for all the fields. This results in an increase of the measured temperatures by 0.9 CKl 


7. REGISTRATION OF NIGHT THERMAL DATA OVER DAY THERMAL AND VISI8LE DATA 


Day thermal and visible data are automatically registered for they were acquired during the same flight through the 
same instrument optics and with the same resolution, 2.5 m at 1000 m altitude. Night thermal data acquired during a pre- 
vious flight, have the same resolution but they must be superposed by software to day data. 


7.1 The method used 

The problem is to nvp the grid of night data on the grid of day data. This is obtained by using bivariate poiy^om a. 's 
as mapping functions [7|. A pixel with coordinates (x, y), given by the line and column position, in the dav data, s thus 
the transformed element from a pixel with coordinates (u, v) in the night data such that: 


u ' u(x, y) 


v - v(x, y) 


* £ V 


© c 0 11*0 


„ £ V 

b^.x ^ 3 

p*0 Q -O 
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where u and v are calculated from x and y. 

In the present case second order polynomials were chosen. The method uses ground control po : nts which are com Ton 
to day and night data to determine the a r , and b. , coefficients by le.ist square fitting. The fitting o' tr.e ground control 
points is performed in the untransformed grid (night data), rather than in the transformed grid (dav data), as wou d seem 
more natural. The latter so’ution would indeed rcc;uire to express the (x, y) coordinates as a function of the (u, v, t.e. to 
reverse the preceding expressions. The calculated fu, v) coordinates will not be generally integer and the nearest ne ghbour 
rule is applied to determine the pixel radiometric value. This simple procedure has been preferred in the present case to 
interpolation or convolution methods, in consideration of the average accuracy of the registration which in any case can 
not guarantee less than two pixels. 

As the complete strip of data is rather tong (more than 4000 lines of 672 pixels), the registration is made for sucoes 
sive blocks of 672 x 300 points with 50 ' overlapping; 12 to 20 control points by block determine the respective oair of 
transformation polynomials. Discontinuities between blocks are avoided by combining with variable weights the coeffi- 
cients cf the polynomials referring to two successive blocks. As the line of data processed moves from the centre of the 
first block to the centre of the second block, the weights relative to the firs pair of polynomials are changed l>e by i ne 
from 1 to 0. the reverse bemg done for the second pair. The process is then repeated with the second and third blocks 
and so on. 




7*2 TlnmuHsobttioad 

They are displayed in Fig. 4. The average registration error between ground control points was found to be 2.3 pixels 
with minimum 0 and maximum 7 pixels. The average registration error on the other points is not known but is very likely 
to have a similar value in this case of low order transformation polynomials. 

a RESULTS - SOIL MOISTURE AND EVAPORATION MAPS 

8.1 Simulated night temperatures 

As already said* the process set up here uses a series of look- up tables generated by the Tel I -us model. For each table, 
Tell-us simulates Td and Tn values which must fit with the various input parameters (weather data, albedo, etc.) a~d with 
a set o< values of THI and RSH. In order to save computing time in the calculation of the 57 necessary tab es, a fast ver- 
sion of Tell-us was run which uses a somewhat high heat capacity to allow calculation time steps of one hour (instead of 
2 to S minutes for the slow version which requires 5 minutes machine time per tablet. As a result, the simu ated T rs va- 
lues are systematically lower by about 1 .5' . A second systematic effect derives from the very low wind speed during the 
night; earlier experiments have indeed shown that the simulated Tn falls down below measured valuas from abo^jj 1 ' in 
this situation. The conclusion was then to increase the simulated Tn values by 2.5 s while waiting for a solut on of :h# 

above problems by improve menu on the calculation scheme and on the model. 

/ 

8.2 Thermal soil moisture 

All mapped results arc relevant to the fields for which Tell*us is presumably applicable. Compare for V jstration 
Figs. 3 to 5 and Figs. 6 to 10. Vegetated zones like meadows, forests, standing cereal crops were not consicered. 

Fig. 9 is the display of both thermal inertia and soil moisture mapped in 6 classes of thermal inertia so i* moisture. The 
soil moisture measurements available for our purpose were done only in the ploughed field and gave a mean value of 0.30 
(.* o ~ 0.024) water content by weight determined over 8 core samples taken at 15 cm depth. The mean va ue of so I 
moisture calculated by Tell-us over the whole field is 0.22 - 0.03 by volume which corresponds to 0.21 r 0.02 by weignt, 
assuming a density of 1,1 for the 15 cm top soil layer* The calculated value is 30^ lower but this comparison, although it 
is indicative for a tendency, should be considered with some caution; first because the moisture content of the core 
sample at 15 cm depth is, not likely to be identical to the moisture content of the 15 cm top layer, second because the 
correspondence between thermal inertia and soil moisture leaves somewhat undetermined the depth of the top iayer con- 
cerned and third because the top layer density is not known exactly and could vary from 1.0 to 1.2. 

A comparison was made with apparent thermal inertia P calculated with NASA-HCMM formula [81. expressed as 
P- CO - a)/(Td - Tn), where a is the albedo and C is a constant which, in the present case, has been determined in 
order to fit with the mean values for the ploughed field of Td, Tn and THI calculated by Tell-us. The results are c-soia' eu 
in Fig. 8 with class limits identical to those used in Fig. 9. As expected, the mapping inside the ploughed f e>d is almost 
identical in the two figures. One can see that the general trends are roughly the same in the other fields. The detai ed va 
nations differ noticeably however, although to a minor degree for the fields towards the left side. Further investigations 
should consequently examine the possibility to approximate the Tell-us model with the help of such a simp e formulation 
under well controlled and well calibrated conditions. 


8.3 Cumulative daily evaporation 

The mapped results obtained with Tell-us are displayed in Fig. 10. It is seen that CDE varies roughly in the same di- 
rection as THI and soil moisture which is reasonable. The detailed behaviour of CDE is however more comp ex because of 
the role played by RSH in the look-up tables. Two fields have for instance higher mean values of CDE than :he three 
fields located between them (towards the left side of Fig. 10), although having about the S3rr,e average THI value as the 
latter; this is due to much higher RSH values, on the average, for the two fields. 

Evaporation is, on the other hand, estimated to be zero (or negative which means condensation) in some zones, for 
8^ of the total area. As condensation did not actually occui during the flight experiment w must conclude that low 
range evaporation is not v/ell calculated by the model. 

Limiting values of CDE have been calculated with the "Tergra" model (91 for volurru a:ure content of 0.44 (i.e. 

THI * 1940) and RSH « 1 .0; they are compared In Table 2 to the Tell-us values which are lower by about 10S. 
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Atwitt be by comp*****? fi^. 9 and 10 witf* Fig. 3. the detailed mapped Hk mm fields *r* mflutfKtd 

Inf the pixel by pixel variation of Td. which is not likely to occur really, as soil moisture arid evaporation will vary more 
slowly. This effect is stronger for non -homogeneous field surfaces trade of variable mixtures of burnt and unburnt stubbles 
or cropped cereals. In these cases there is also a strong variation of aerodynamic roughness on a pixel by pixel basis. The 
model is very sensitive to these vena Jons: doubling the roughness has the effect to modify by 1 to 3' the simulated day 
and night temperatures. In such, difficult situations, the maiming of moisture and evaporation would then require that 
further to attaedo. also roughness be acquired and taken into account on a pixel by pixel basis, which is hardly fees bfe at 
the present stage. A possible solution could be to smooth the effect by averaging over superpixels, if one agrees, ho .sever, 
to the corresponding loss of detailed information. 


9. CONCLUSIONS 

From the point of view of data processing it has been verified that a detailed mapping of thermal inertia, soil moisture 
and evaporation is possible with day and night thermal and albedo data acquired by air borne scanner, and by using a 
rather complex model like "Teil us" which reverses a temperature simulation process. 

The comparison with measured values is subject to uncertainties but it seems to indicate a tendency of the model to 
underestimate moisture content. Other measurements would be necessary to check the results on a more re : able basis. 

The same can be said for cumulative daily evaporation for which no measurement was available; the above mentioned 
tendency *s however afd verified. 

In inhomogeneous areas, where temperature, albedo and surface roughness vary rapidly, these rapid variations tend 
to mask partially the more smoothed evolution of moisture and evaporation, as long as roughness is not taken into 
account pixel by pixel, as was done for albedo. 
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FIG.1* 

Look-up graphs calculated by in the cm* o* bare sod with albedo 0.09 


Correction eu'vet :o rerro%e :hr 
scan angle effect cn loecvc-netrc 
vatucv Upper curvet: vr>r. m c' 
the coiumn mean and r*’a: ve fitted 
turve along the Kan i-ne. 7n* u*e 
jr variances in the lev**' c.'vn. 

I our channeitonly a^t rec-esroted: 
day thermal (2a). mgnt mp'rrval (2bl 
channel 5 (2c) and channel 10 (2o). 
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Fig. 4a: Mapped night temperatures, ori- 
ginal data. Range of variation 
O'C 18"C (radiative values). 


Fig. 4b: Mapped night temperature*. 

Data geometrically corrected. 
The transformation is visible 
through the grid deformation 


Fig 3 : Mapped day temperature Range of 
of variation 8°C - 44°C (radiative 
values) 


Fig. 5 Mapped albedo values Range of 
variation 0 02 0.20 in the in- 

vestigated fields 
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0 0.14 0.27 0.40 0.53 

F:g 10 Cumulative daily evapora 
lion calculated bp Tel! us 
and mapped n 5 grey 
levels (mm) 1: bare soil: 

2 stubbles. 3 burnt & 
unburnt stubbles 4: stub 
bles w»th green shoots 
5: stubbles & cropped 
cereals. 

Thermal inertia 
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Fig. 9. Thermal inertia and soil 
moisturp calculated by 
Tell-us and mapped in 6 
grey levels. 


900 1150 1400 1650 1900 2150 2400 


Fig. 8 Apparent thermal inertia 
mapped in 6 grey levels 


Fig G Temperature diffeiences 

(Td - Tn) mapped in 6 grey 
levels 
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Fig. 7 DiffLMence between Fig. 9 and 
Fig. 8 mapped in 5 grey levels. 





